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ABSTRACT 

With the high-resolution spectroscopy now available in the optical and satellite UV, it is possible to determine 
the neutral/ionized column density ratios for several different elements in a single cloud. Assuming ionization 
equilibrium for each element, one can make several independent determinations of the electron density. For 
the clouds for which such an analysis has been carried out, these different estimates disagree by large factors, 
suggesting that some process (or processes) besides photoionization and radiative recombination might play an 
■ important role in the ionization balance. One candidate process is collisions of ions with dust grains. 

Making use of recent work quantifying the abundances of polycyclic aromatic hydrocarbon molecules and other 
grains in the interstellar medium, as well as recent models for grain charging, we estimate the grain-assisted ion 
£> ' recombination rates for several astrophysically important elements. We find that these rates are comparable to 

the rates for radiative recombination for conditions typical of the cold neutral medium. Including grain-assisted 
ion recombination in the ionization equilibrium analysis leads to increased consistency in the various electron 
density estimates for the gas along the line of sight to 23 Orionis. However, not all of the discrepancies can be 
eliminated in this way; we speculate on some other processes that might play a role. We also note that grain- 
assisted recombination of H + and He + leads to significantly lower electron fractions than usually assumed for the 
cold neutral medium. 

Subject headings: dust — ISM: abundances — ISM: clouds — stars: individual (23 Orionis, HD 215733) 
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I | 1. INTRODUCTION Lepp et al. 1988). The "removed" ion might remain in the gas 

i-^ , ™ , . , , .. ,- cc ■ . . , , phase (in a lower ionization stage) or it might be depleted from 

*H The electron number density n, in diffuse interstellar clouds \ , , . , b , . f _ r . . 

> 1 . £ , ■ c , c .1 * j i , c a cc the gas altogether, if it sticks to the grain. The grain-assisted 

c/2 , is often inferred from the observed column densities of differ- b f '. b r P 

rr* , ■ ■ u . c . , T , .■ c ■ ■ .. ion removal rate is expected to vary from species to species 

ent ionization stages of metals. The assumption of ionization . , / . / , \ / 

• • . • ■ , , since the possibility of charge transfer depends on the electron 

eouilibnum vields 

. £L ' affinity of the ion and the collision rate depends on its mass. 

r(X')n(X') = a r (X' +1 , r)n c n(X' +1 ) , (1) The collision rate also depends on the grain charge and size; 

j_J ■ when the grain charge Z < 0, the ratio of the collision cross sec- 

£ where »(X<) is the number density of the i-th ion of element X, ^ tQ ^ tric cross section is l for smaller ins 

T(X') is the ionization rate for X' X' +1 , and a r (X' +1 , T) is the (Draine & Sutin 198?) Jhu ^ a j lation of moleC ules 

rate coefficient (which depends on the gas temperature T) for Qr smaU ing wkh z < Q CQuM contribute significant i y 

recombinations X — * X . It is usually assumed that radiative t j Qn remova j 

recombination is dominant at low temperatures; dielectronic re- Observations of emission from the diffuse interstellar 

combination is often important when T > 10 K. medium (ISM) haye revealed features at 3 j fi 2> ? n %£ and 

The Goddai-d High Resolution Spectrograph on the Hubble 1 1 3 (see SeU 1994 for a reyiewX which haye been iden . 

Space Telescope has made possible the simultaneous determi- tified as C _ H and c _ c stretchin and bendin modes in j 

nation of « f through analysis of the ; ionization balance of several Hc hydrocarbons (PAHs ; Leger & Puget 1984; Al- 

different metals. In some cases, the values of n e inferred from lamandola> TielenSi & Barker 1985) . Li & Draine (2001) com . 

different species differ by large factors (Fitzpatnck & Spitzer red observations of diffuse Galactic emission with detailed 

1997; Welty et al. 1999), suggesting either that the ionization is model calculations for PAHs heated b Galactic starli ht and 

far from a steady-state equilibrium or that another process (or found ^ a c abundance 2 _ 4 _ 6 x 1Q -5 is r ired in hydrocar . 

processes) besides photoionization and radiative recombination bon molecules with < 1Q 3 c atoms .3 In recent charging modelSi 

affects the ionization balance. PAHs are found to be predominantly neutral in the diffuse ISM 

One candidate process is ion recombination via electron (Bakefj & Tielenf , 1994; Wei tner & Draine 2001b) . 

transfer from a grain or molecule to the ion during a col ision In R wg estimate rain . assisted ion removal rates for sev . 

(Weisheit & Upham 1978; Omont 1986; Draine & Sutin 1987; 5 

1 Henceforth, we will include molecules in the "grain" population. 

2 By "abundance" we mean the number of atoms of an element per interstellar H nucleus. 

3 Electric dipole emission from this population can account for the dust-correlated component of the diffuse Galactic microwave emission (Draine & Lazarian 1998; 
Draine & Li 2001). 
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eral metals, employing a grain size distribution that (a) includes 
PAHs in numbers sufficient to account for the observed infrared 
and microwave emission and (b) reproduces the observed ex- 
tinction of starlight (Weingartner & Draine 2001a). In §3, we 
compare the grain-assisted ion removal rate with the radiative 
recombination rate, finding them to be comparable in the cold 
neutral medium. In §4, we discuss the determination of elec- 
tron densities from observed column density ratios when grain- 
assisted ion removal is important. In §5, we apply these re- 
sults to two observed lines of sight. We find that the inclusion 
of grain-assisted ion removal cannot, on its own, reconcile the 
electron densities inferred from different species observed to- 
wards 23 Ori. In §6, we speculate on additional processes that 
might affect ionization balance along this line of sight. In §7, 
we investigate the impact of grain-assisted ion recombination 
on the electron fraction in the cold neutral medium. We sum- 
marize our conclusions in §8. 

2. ION REMOVAL RATES 

Consider the collision of an ion X' with a dust grain of ra- 
dius a and charge Ze (e is the proton charge). We assume 
that an electron is transferred from the grain to the ion if 
IP(a,Z) < IP(X' _1 ), where IP(a,Z) is the grain ionization po- 
tential and IP(X' _1 ) is the ionization potential of X' -1 (equal to 
the electron affinity of X ! ; see Table 1 for our adopted values). 
This is certainly an oversimplification since the electrostatic en- 
ergy of the grain-ion combination usually increases following 
electron transfer. However, it is not clear how to evaluate the 
electrostatic energy when the distance from the grain to the ion 
is < a few angstroms. Thus, our estimates for ion removal rates 
should be regarded as likely upper limits. Below, we will es- 
timate the extent to which these rates might be changed if the 
electrostatic energy changes were included in the analysis. 

We adopt the following expression for the ion removal rate: 

1 dn^ffaZViWeUPQC-^-IPiaM] , (2) 

z 



dt 



where n, is the number density of X ! , n gr (a) is the number den- 
sity of grains with radii < a, f(a,Z) is the fraction of the grains 
of radius a that have charge Ze, and 9(y) = if y < and 1 if 
y > 0. Ions arrive at a grain with size a and charge Z at the rate 



Ji(Z) = m 



8kT 



Km, 



1/2 



A^ 2 Ka z J(Ti = akT/qf,£i = Ze/qf) 



(3) 



where k is Boltzmann's constant, m p is the proton mass, A x is 
the ion mass number, qi is the ion charge, and expressions for 
the dimensionless function /(t;, £,-) are given in Draine & Sutin 
(1987). 

Grain charge distributions in the diffuse ISM are set pri- 
marily by the balance between electron loss via photoelectric 
emission versus electron gain via accretion from the gas. We 
calculate the steady-state charge distributions f(a,Z) using the 
charging algorithm given by Weingartner & Draine (2001b) 
and the average interstellar radiation field (ISRF) spectrum for 
the solar neighborhood as estimated by Mezger, Mathis, & 
Panagia (1982) and Mathis, Mezger, & Panagia (1983). An 
often-used measure of the radiation intensity is the parameter 
G = «rad/ M Hab' wriere «"ad is the energy density in the radiation 
field between 6 eV andl3.6eV and M^ v ab = 5.33 x 10" 14 ergcm" 3 
is the Habing (1968) estimate of the starlight energy density in 

4 For comparison, the interstellar radiation field estimated by Draine (1978) has u = 



this range. 4 For the ISRF, G = 1.13. Grain charging is largely 
determined by the parameter 

il> = GVf/n e (4) 
(Bakes & Tielens 1994; Weingartner & Draine 2001b), with 
only a weak additional dependence on T. Examples of charge 
distribution functions f(a,Z) are given in Weingartner & Draine 
(2001b). 

We define a rate coefficient a g (X',ip, T) for electron transfer 
to ion X ! due to collisions with grains by 



drii • 

— =-a g (X ,tp,T)nin H 



(5) 



where «h is the H nucleus number density. In Figure 1, we dis- 
play a g as a function of ip for several ions X' and T = 100K. In 
Figure 1, and throughout this paper, we employ the grain size 
distribution from Weingartner & Draine (2001a) with R v = 3.1 
and be = 6 x 10" 5 (their favored distribution for the average dif- 
fuse ISM). Note that a g decreases sharply as ip increases, i.e., 
as the average grain charge increases. When T < 10 3 K, the de- 
pendence of a g on T is quite mild since, for neutral grains and 
t, <C 1, J oc r -1 / 2 , cancelling the factor of T l l 2 in equation (3). 

For the cases illustrated in Figure 1, we have also computed 
a modified ion removal rate coefficient a' g , for which we apply 
a different criterion for electron transfer: 

IP(a,Z)-IP(X'- 1 ) + Af/(Z,0<0 . (6) 
Here, AU is the change in the interaction energy between the 
grain and the ion due to the electron transfer. We approximate 
(see Appendix) 

A^(Z,z) Mfe -Z-l)-f- + - **- l £f + - , (7) 
(a + ro) 2ro(fl + ro) z (2a + r ) 

where Ze is the initial grain charge, ze is the initial ion charge, 

and ro is the atomic radius (see Table 1 for our adopted values 

of ro). Since the distance of closest approach between the ion 

and the grain surface is unknown, we have taken it to be equal 

to the atomic radius. 

The first term in equation (7) is the monopole contribution 
due to the net charges on the grain and ion. Since the grains 
are predominantly neutral (Z = 0) and most ions of interest are 
singly ionized (z = 1), the inclusion of this term usually has a 
negligible effect on the calculated ion removal rate. The second 
term in equation (7) is due to the polarization of the grain by 
the ion. Our classical calculation of the interaction energy as- 
sumes that the grain is (a) a perfect conductor and (b) a sphere 
with a perfectly sharp surface. These assumptions can lead to 
substantial errors in the estimate of the image potential when 
the ion-grain surface separation < a few angstroms. We have 
neglected the dipole moment in the atom or ion induced by the 
grain charge, since the approximation of the atom/ion dipole as 
a point dipole breaks down when the distance from the atom/ion 
to the grain surface is comparable to the size of the atom — the 
case of interest here. 

In Figure 2, we display a' g /a g . For cold and warm neutral 
medium conditions, a' g /a g is generally <~ 1; thus, it appears 
that, for most cases of interest, an accurate accounting of the in- 
crease in electrostatic energy following electron transfer would 
not substantially reduce the computed ion removal rate. 

For a given ion, a g can be represented approximately by the 
following fitting formula: 

nn I t\ 10- 14 C cm 3 s- 1 

ag( } ~ l+QV^ [i+c 3 T*iirCr<** T ] ' ( j 

8.93 X 1(T 14 ergenr 3 between 6 and 13.6 eV, or G = 1.68. 
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FIG. 1 . — Rate coefficient a g for grain-assisted removal of selected ions as a function of the charging parameter ip, for the ISRF spectrum and gas temperature 
r = 100K. 
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FIG. 2. — a' g /a g , where a' g is computed assuming a more stringent criterion for the transfer of an electron from grain to ion, in which the increase in electrostatic 
energy due to the interaction of the grain and ion charge is taken into account in an approximate manner. Nominal cold and warm neutral medium conditions are 
indicated ("CNM" and "WNM", respectively). 



with T in kelvins and ip in units of K 1,/2 cm 3 ; values of the 
fitting coefficients Co — C(, are given in Table 2 for selected 
ions. 5 Approximation (8) is accurate to within 20% when 
10 < T < 10 3 K and 10 2 < ip < 10 6 K 1 / 2 C m 3 . In the grain- 
assisted recombination of Ca ++ , there is a possibility for two 
electrons to be transferred to the ion, resulting in Ca°. We as- 
sume that this occurs in a two-step process and that the energy 
released following the first electron transfer is dissipated into 

5 A FORTRAN subroutine implementing approximation (8) is available on the 



heat, and thus is not available to assist in the second transfer. 
With these assumptions, the criterion for double electron trans- 
fer is 

IP(a,Z+l)<IP(Ca°) . (9) 

For this criterion, 10 < T < 10 3 K, and 10 2 < tp < 10 6 K 1 / 2 C m 3 , 
the probability of double charge transfer can be approximated 

Wide Web at www.cita.utoronto.ca/~weingart. 



4 



Weingartner & Draine 



to within 10% by 

fcs 



0.415 



*Ca° : 



1+6.38 x 10- 5 ?/> L252 (l + 3.57 x lO" 6 ^ 1124 ) 



(10) 

Note, however, that the double electron transfer rate could 
be dramatically reduced if electrostatic energy changes (as in 
eq. [7]) are included in the criterion for electron transfer. 

3. GRAIN-ASSISTED ION REMOVAL VERSUS RADIATIVE 
RECOMBINATION 

To gauge the importance of grain-assisted ion removal in the 
diffuse ISM, we plot in Figures 3 and 4 a critical electron frac- 
tion 

a g QC,il>,T) 



Xait (X l ,lP,T): 



a r (X',T) 



(11) 



if x = rie/rin < x cl n, then grain-assisted removal of X' occurs 
more rapidly than radiative recombination to X'" 1 . We use the 
FORTRAN routine rrfit.f to evaluate radiative recombination 
coefficients. 6 In Figure 3, we take T = 100 K, appropriate for 
the cold neutral medium (CNM). For a nominal CNM electron 
fraction x w 10~ 3 and tp w 400K 1/2 cm 3 , the grain-assisted ion 
removal rate is generally comparable to the radiative recombi- 
nation rate and therefore will appreciably affect the ionization 
balance for H + and most metal ions. 

In Figure 4, we take T = 6000 K, appropriate for the warm 
neutral medium (WNM) and warm ionized medium (WIM). 
Adopting x w 0. 1 and ip w 3000K 1 / 2 cm 3 for the WNM, we find 
that the grain-assisted ion removal rate < 10% the radiative re- 
combination rate. For the WIM (x«l,^« 1000K 1/2 cm 3 ), 
grain-assisted ion removal contributes at most at the few per- 
cent level. 

4. ELECTRON DENSITY DETERMINATIONS INCLUDING 
GRAIN-ASSISTED ION REMOVAL 

When grain-assisted ion removal is important, the ionization 
balance equation (1) must be replaced with the following set of 
equations: 

dn M 



r(X')«,- a r (X i+l )n e ti i+1 - aJX' +l )n H n i+l 



dt 



(12) 



and 

^=-T(X i )n i + a r (X i+l )n e n i+l + a g (X i+l )n li n i+l (l-s) , (13) 

where s is the probability that an ion sticks to a grain following 
charge exchange. 7 Of course, if more than two ionization stages 
are populated, then more equations are needed. When s > 0, n,- 
and both decay, but, with the exception of transients that 
die out rapidly, the ratio is constant: 



2b 



1+1 + 



4b 



1/2 



(14) 



where b = r\ + r 2 (l — s), q = \ — r\ - r 2 , r\ = a r n e /Y, and 
r 2 = cym/r. 

Solving equation (14) for n e in terms of the observed column 
density ratio R = n(X i )/n(X i+1 ), we find 

RT 



n e 



1- 



l—s+R a g tin 
l+R RT~ 



(15) 



Since a g depends on n e (through its dependence on the charging 
parameter ip), equation (15) must be solved iteratively, but gen- 
erally converges rapidly. Usually ri<l and r 2 <C 1 (implying 
R <C 1). Thus, the impact of grain-assisted ion removal on the 
electron density determination decreases as the sticking proba- 
bility s — ► 1 . The observed interstellar depletions vary widely 
from element to element, implying s< 1 for some elements 
(e.g. S) while s ~ 1 for others (e.g. Ti, Ca; see Weingartner 
& Draine 1999). This wide variation in s could also lead to 
discrepancies in electron density estimates when different ele- 
ments are observed. 

A large fraction of the Ca in diffuse clouds can be in Ca ++ ; 
thus, we must consider the ionization equilibrium of a 3- 
level system in this case. For arbitrary s, this results in a 
rather complicated cubic equation for n e in terms of the ob- 
served ratio R\ = n(Ca°)/n(Ca + ). However, the equation sim- 
plifies to a quadratic when s = or s = 1. Before display- 
ing the results, it is convenient to introduce some notation: 
n e = /?ir(Ca°)/av(Ca + ) (this is the value of n e when grains are 
unimportant), g\ = a i ,(Ca + )/a r (Ca + ), g 2 = a g (Ca ++ )/o! r (Ca ++ ), 
d = ov(Ca ++ )/a r (Ca + ), 7 = r(Ca + )/r(Ca°), h = 
a g (Ca + )n H /a r (Ca + )n e , h 2 = a g (Ca ++ )n H / a r (Ca + )n e , and f 2 = 
/ca^^ca" (fromeq. 10). 

When 5 = 0, 



1 1 

•= 2(^"^i"h)(1 + Go)-2^2«h(1-2o) 



where 



Qo=U 



^■ginHjfifie 



R\[(.gi-g2)nn-h e ] 2 



1/2 



(16) 



(17) 



Note that equation (16) reduces to equation (15) when 7 = or 
f 2 = 0. The ratio R 2 = n(Ca°)/«(Ca ++ ) is given by 



R 2 = 



h 2 +d n 4 

n, 



(18) 



When s=l, 



n e = n e 



where 



RxhiQ. +dR x ){\ + Q x ) + R l {d+h 2 ){\ +J?0(1 -20 + 7 



2{\+Ri){\+dR{) 



(19) 



Gi 



1 + 



j 2 + 2jR l {h l (l+dR l ) + (l+Ri)[d-h 2 (l+2dR l )]} 



R\[h l (\+dR{)-(d+h 2 ){\+R l )} 2 



1/2 



Ri = 



Ri 

1 



h 2 -R~ 1 1 + (d+R~ 1 



(20) 
(21) 



6 The subroutine rrfit.f was written by D. A. Verner and is available on the World Wide Web at www.pa.uky.edu/ 

7 We assume that ions do not stick to grains following collisions without charge exchange. 



Note that equation (19) reduces to equation (15) when 7 = 0. 

5. APPLICATION TO OBSERVATIONS 

5.1. HD 215733 

Fitzpatrick & Spitzer (1997) performed ionization equilib- 
rium analyses of Ca + /Ca ++ , C°/C + , Mg°/Mg + , and S°/S + for 
several clouds along the line of sight to HD 215733. In the 
case of Ca + , the Ca" 1 " 1 " abundance was estimated based on an 
assumed gas-phase Ca abundance. In Table 3, we reproduce 
their inferred temperatures and H number densities for the cold 

verner/fortran.html. 



Recombination on Grains 



5 



10- 



10- 3 



10 4 - 



10- 5 



10- 6 




T=100 K 



ISRF spectrum 



Na± \_\ \ 



CNM 



Ca 



10 10 2 10 3 10 4 10 5 

cm- 1 ) 



10 6 



FIG. 3. — Critical electron fraction: if x < x a n, then the indicated ion is removed more effectively via collisions with grains than via radiative recombination. The 
ticks labeled "CNM" are for canonical cold neutral medium values of ■)/> and x. 
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FIG. 4. — Same as Figure 3, except for T = 6000K instead of T = 100 K. The ticks labeled "WNM" are for canonical warm neutral medium values of ip and x. 



cloud components. In the upper panel of Figure 5, we repro- 
duce their electron density determinations. 8 The electron den- 
sities inferred from the C and Mg abundance ratios appear to be 
roughly consistent, but those inferred from the Ca abundance 
ratios are systematically lower. 

In the lower panel of Figure 5, we display revised n e deter- 
minations that take grain-assisted ion removal into account. We 
assume s«l for Ca, which is observed to be highly depleted 



in cold clouds, and skO for C, Mg, and S, which are not as 
severely depleted. With this modification, the Mg determina- 
tions lie between the C and Ca determinations, and are consis- 
tent with both. The C and Ca determinations remain inconsis- 
tent. If dissociative recombination of CH + is important in these 
clouds, then the n e inferred from the C abundance ratios are 
too high; unfortunately, the CH + column density towards HD 
215733 is not known. Alternatively, the n e inferred from the 



8 The Mg + f-values used by Fitzpatrick & Spitzer (1997) are apparently too large by a factor of Ri 2.4 (Fitzpatrick 1997). Thus, we have multiplied the n e values 
obtained from the Mg column density ratio by 2.4. 
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Component Number 

FIG. 5. — Upper panel: Fitzpatrick & Spitzer (1997) results for n e (in cirT 3 ) for cool cloud components along the line of sight to HD 215733, with n e deter- 
mined from ionization equilibrium for Ca + , C°, Mg°, and S° (Mg results are adjusted as per Fitzpatrick 1997). Lower panel: Modified n e determinations, with 
grain-assisted ion removal included. We assume s = 1 for Ca and s = for C, Mg, and S. 



Ca abundance ratios could be too small, as would be the case if 
the Ca abundance has been overestimated. Unfortunately, with 
only three elements for which both ionization stages have been 
observed, it is difficult to arrive at definitive conclusions regard- 
ing the electron densities in many of these components. 

5.2. 23 0ri 

Welty et al. (1999, hereafter W99) have conducted a de- 
tailed investigation of the diffuse clouds along the line of sight 
to 23 Orionis. They inferred T ~ 100K, «h ~ 10cm" 3 , and 
Nh ~ 5 x 10 20 cm" 2 for the cold cloud material along this sight- 
line at heliocentric velocities 20 < v Q < 27kms" 1 . They also 
inferred n e using equation (1) for ionization equilibrium (be- 
tween X° and X + ) with several different elements, and found 
widely varying values; we display their electron densities as 
open triangles in Figure 6. 9 The electron densities inferred if 
ionization equilibrium is assumed vary over about an order of 
magnitude, depending on which element is used. 

Electrons will collisionally excite the 2 P3/2 excited fine struc- 
ture levels of C + and Si + . W99 report N(C + ) = 10 16 95 cm" 2 , and 
N(C + *) < 10 15 cm" 2 , from which we infer n e < 0.075 cm" 3 ; we 
assumed T = 100 K, a collision strength O = 2.1 from Keenan 
et al. (1986) and H collisional rates from Launay & Roueff 
(1977). 

For Si\ W99 find N(Si + ) = 10 15 ' 37 cnr 2 and N(Si + *) < 
10 la7 cm" 2 , from which we find n e < 0.12cm" 3 , assuming a 
collision strength £1 = 5.74 (Keenan et al. 1985). H collisional 
rates have not been calculated; we use the rate coefficient cal- 
culated by Launay & Roueff for deexcitation of C + ( 2 P 3 / 2 ). 

Upper limits on n e derived from fine-structure excitation are 
shown in Figure 6. 

The electron densities inferred from ionization equilibrium 

9 We display the W99 n e determinations using photoionization rates for the "WJT 
for the ISRF of Mathis et al. (1983). 



including grain-assisted ion removal are displayed as filled tri- 
angles in Figure 6. We have assumed s = in order to maximize 
the effect that grain-assisted ion removal might have. (For most 
of the elements in Figure 6, this assumption seems reasonable. 
For Fe and Ca, however, the large depletions observed in inter- 
stellar gas suggest that s > 0.1.) The W99 electron densities in- 
ferred using C, Na, Mg, and K are particularly large; the result- 
ing n e j riu are high enough that the inclusion of grain-assisted 
ion removal does not dramatically alter the inferred values of 
n e . Thus, the various n e determinations do not appear to be 
reconciled by including this process. 

However, some of the W99 n e determinations might be too 
high for other reasons. Dissociative recombination of CH + 
might be an important additional recombination mechanism for 
C. W99 found that if the CH + observed towards 23 Ori re- 
sides in the cold cloud for which the C° and C + column den- 
sities were measured, then the estimated n e must be revised 
downward from 0.26 to 0.04cm" 3 (open square in Figure 6); 
of course, grain-assisted ion removal would reduce it further 
(filled square). 

Na + and K + are unobservable since they have no resonance 
lines below the Lyman limit. Thus, W99 had to assume gas- 
phase Na and K abundances in order to analyze the ionization 
balance for these elements. W99 assumed depletions of 0.5 dex 
(relative to solar) for both elements. If, instead, these elements 
are assumed to be undepleted, then the n e estimates are substan- 
tially smaller; these are indicated in Figure 6 by filled (open) 
squares for analyses that do (do not) include grain-assisted ion 
removal. 

Note that Welty & Hobbs (2001) argue, on the basis of a large 
number of observed sight lines, that Na and K are generally de- 
pleted from the gas phase by 0.6-0.7 dex relative to solar. We 
consider the depletion estimates for Na and K to be quite uncer- 

radiation field (see Table 1); these photoionization rates are very similar to those 
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FIG. 6. — Electron density n e inferred from ionization equilibrium for the cold gas towards 23 Orionis. Open triangles: Welty et al. (1999) analysis; Filled trian- 
gles: Welty et al. (1999) analysis plus grain-assisted ion removal; Open squares: Modified analysis (see text), without grain-assisted ion removal; Filled squares: 
Modified analysis, including grain-assisted ion removal. For each species, our favored n e determination is indicated by an enlarged symbol. Our actual favored 
n e 0.03 cm"" 3 is indicated by the dashed line. Upper limits from the observed fine structure excitation from C + and Si + are also indicated. 



tain, as they rely on knowledge of the electron density (which, 
as we have seen, is uncertain since "standard" ionization anal- 
yses give discrepant results when multiple elements are used). 

If the revisions described above (for C, Na, and K) are cor- 
rect, then the observed column densities for six elements (C, 
Na, S, K, Mn, Fe) yield n e k, 0.03 cm" 3 , consistent with the 
upper limits n e < 0.075 cm" 3 and n e < 0. 12cm" 3 obtained from 
the fine structure excitation of C + and Si + . For « c /«h = 3 x 10~ 3 , 
grain-assisted ion removal contributes significantly to the re- 
combination. 

However, discrepancies remain for three elements: Ca, Mg, 
and Si. 

1 . W99 found n e w 0.95 cm" 3 from Ca°/Ca + ionization bal- 
ance. Using equation (16), we find n e = 0.93cm" 3 even 
if s = 0, so including recombination on grains does not 
resolve this discrepancy. 

2. Even after allowing for grain recombination with s = 0, 
the Mg°/Mg + abundance ratio implies an electron den- 
sity n e rts 0.2cm" 3 — a factor of 6 greater than our best 
estimate. 

3. Even after allowing for grain recombination with 5 = 0, 
the Si°/Si + abundance ratio implies an electron density 
n e w 0.08 cm" 3 — a factor of 2.5 greater than our best es- 
timate. 

Unless the photoionization rates have for some reason been 
overestimated, the observations suggest that additional physical 
processes are involved for these three elements. 

6. POSSIBLE ADDITIONAL PROCESSES AFFECTING 
IONIZATION BALANCE IN 23 ORI 

6.1. Dust Destruction ? 



In Figure 7, we plot the inferred electron density, normalized 
to our favored value of 0.03 cm" 3 , versus the gas-phase element 
abundance in the cold cloud towards 23 Ori, normalized to the 
standard cold cloud abundance (see Table 5 in W99). We have 
included two points for Fe, with sticking coefficient s = and 
s = 1 ; Fe is generally strongly depleted in cold clouds, suggest- 
ing that the sticking coefficient may be ~ 1. Figure 7 shows 
a clear trend: greater neutral abundance enhancements are ob- 
served for elements with greater overall abundance enhance- 
ments. The total gas-phase abundances of Na and K are not 
known, but if they are enhanced by w 0.4 dex over their stan- 
dard values, then they fit the trend nicely. (We have also in- 
cluded points for Na and K for which the gas-phase abundance 
is enhanced by 0. 1 dex, as assumed by W99, and by 0.6 dex, as 
assumed by us above. Neither of these choices lie close to the 
trend.) The observed correlation suggests the possibility that 
the enhanced neutral Mg, Si, and Ca abundances result from 
ongoing destruction of dust containing these elements, with the 
atoms injected into the gas as neutrals. However, the production 
rate of neutrals due to dust destruction would have to exceed the 
rate due to recombination; it is hard to imagine how such rapid 
dust destruction could be occurring. It seems to us unlikely that 
dust destruction is affecting the ionization balance observed to- 
ward 23 Ori. 

6.2. Dielectronic Recombination? 

It is intriguing to note that Mg° (3s 2 ground state), Si (3p 2 
ground state), and Ca° (Ap 1 ground state) are all species for 
which dielectronic recombination is important for T < 10 4 K; 
the dielectronic recombination rate for Mg + — > Mg° exceeds the 
radiative recombination rate for T > 5800 K. W99 argue for a 
temperature T w 1 00 K for this velocity component, but perhaps 
there is some gas in this velocity range with 5000 < T < 10 4 K 
and a relatively high electron density, leading to preferential 
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recombination of those species with the largest (radiative + di- 
electronic) rate coefficients at these temperatures. It is also in- 
triguing to note that W99 report detection of Al ++ at the velocity 
of the dominant neutral absorption, suggesting that some of the 
material in this velocity range may be photoionized by pho- 
tons beyond the Lyman limit, with accompanying ionization 
of H and presumably enhanced values of n e and temperatures 
T s=s 10 4 K. Recall that CH + is present in this velocity range. 
The origin of interstellar CH + remains unclear, but all current 
scenarios for CH + formation, whether in MHD shock waves 
(Draine & Katz 1986; Pineau des Forets et al. 1986), turbulent 
boundary layers of clouds (Duley et al 1992), strong Alfven 
waves (Federman et al. 1996), or intermittent dissipation of tur- 
bulence (Falgarone, Pineau des Forets, & Roueff 1995; Joulain 
et al. 1998), involve regions of enhanced temperature. 

We have found that if w 10% of the column density arises 
in a warm component with T w 10 4 K and n e w 4 cm" 3 , then 
the predicted column density ratios for all of the elements in 
Figure 6 are within a factor of w 2 of the observed ratios. In 
this case, however, the predicted populations of the excited fine 
structure levels of C + and Si + substantially exceeds the observa- 
tional upper limits (see §3.3.3 of W99). We were unable to find 
any models for which (a) dielectronic recombination in warm 
gas significantly increases the neutral/ionized column density 
ratios of Mg, Si, and Ca and (b) fine structure levels are not 
overpopulated. We conclude that dielectronic recombination 
cannot explain the observed high values of Mg°/Mg + , Si°/Si + , 
andCa°/Ca + . 10 



ing for the lower n e inferred from these species. W99 noted 
that charge transfer rates for T w 100K have not been calcu- 
lated for these elements with protons, although potential en- 
ergy curves calculated for SH + (Kimura et al. 1997) imply 
that the rate coefficient for S + H + — > S + + H is very small at 
low temperatures, so that protons would not affect the S°/S + 
ionization balance. It should also be noted that with the limit 
n e < 0.075 cm -3 from C + fine structure excitation, the charge 
transfer ionization rates n(H + )(av) <2x 10" 10 s _1 for Mn and 
Fe since (av) <3x 10~ 9 cm 3 s _1 even if MnH + or FeH + have 
favorable potential energy curves. Charge transfer with protons 
is therefore at most comparable to photoionization for Mn and 
Fe (see Table 1), so that n e inferred from Mn and Fe cannot be 
raised by more than a factor of 2 even if charge transfer with 
protons is maximally effective. We therefore believe that the 
electron densities inferred from S, Mn, and Fe are reliable. 

It is intriguing to note that the three atoms with np 2 ground 
states [C(2p 2 ), Si(3/? 2 ), and Ca(4p 2 )] are all overabundant. In 
the case of C we can attribute the overabundance to dissociative 
recombination of CH + , and the question arises whether disso- 
ciative recombination of SiH + and CaH + might also be taking 
place. CH + is believed to be produced by the exchange reaction 
C + + H 2 + .40e V -> CH + + H. The reaction X + + H 2 -» XU + + H 
is much more endothermic for Si + (1.29 eV) than for C + (0.40 
eV). While the scenario for CH + formation remains unclear, the 
increased endothermicity for Si (and presumably for Ca, though 
the heat of formation of CaH + is unavailable) makes such a 
chemical explanation appear unlikely. Furthermore, production 
of SH + would appear more favorable than, say, SiH + , yet the 
observed S°/S + ratio can be understood without invoking disso- 
ciative recombination of SH + . 



6.3. Chemistry? 

W99 suggested that n e w 0.15 cm -3 in the cold gas towards 
23 Ori and that charge exchange with protons might be an 
important ionization mechanism for S, Mn, and Fe, account- 

10 also conclude that dielectronic recombination is not likely to be responsible for the enhanced Ca°/Ca + . They note that the n e derived from Ca°/Ca + are often 
higher than those derived from other neutral/ionized ratios, but that the line widths of Ca° usually indicate an origin in cool gas. 
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6.4. Atomic Physics? 

The overabundance of the np 1 species suggests that some 
common atomic physics might be involved. If, for example, 
photoionization rates have been overestimated, or radiative re- 
combination rates at T w 100K have been underestimated, the 
rates for all three elements (C, Si, Ca) might be similarly af- 
fected. The photoionization rates are probably unlikely to be se- 
riously in error, but possible underestimation of the low temper- 
ature radiative recombination rates should be considered. The 
Milne relation for the rate coefficient for radiative recombina- 
tion X + + e ~ >X+hv is 

A 1 />00 

(22) 

where the sum over j is over all of the terms for the ground con- 
figuration (e.g., 2s 2 2p 23 Y>, 1 D, ] S for C) and all of the terms for 
one-electron excited configurations of the atom (e.g., 2s2p 35 S, 
2s 2 2p3s 3 P, 2s 2 2p3s l P, 2s2p 33 D, ... for C), g x + is the degen- 
eracy of the ion X + , gj is the degeneracy of excited state j of 
the atom X, Ij is the energy required to ionize from excited 
state j, and <j pi j(hv) is the photoionization cross section from 
excited state j. Equation (22) shows that the radiative recom- 
bination rate at T < 100K could be large if some state j has 
a photoionization cross section a P ij with a resonance within 
< 0.01 eV of threshold. While we have no indication that such 
a resonance exists, we encourage atomic physicists to reexam- 
ine the calculated photoionization rates and implied radiative 
recombination rates. Such a near-threshold photoionization res- 
onance, and consequent large low-temperature recombination 
rate, could conceivably also occur for Mg°. 

7. ELECTRON FRACTION IN THE COLD NEUTRAL MEDIUM 

In this section, we investigate the extent to which grain- 
assisted removal of H + and He + decreases the electron fraction 
x e in the CNM. There are two separate contributions to the total 
electron fraction. The first is due to the ionization of metals, 
primarily C, by the far-ultraviolet radiation field. Since we ex- 
pect C to always be primarily in C + in cold diffuse clouds, we 
simply set this contribution to x e p ~ 2 x 10~ 4 . The second con- 
tribution is due to ionization (primarily of H and He) by cosmic 
rays and by the diffuse X-ray radiation field. 

Wolfire et al. (1995) estimate the cosmic ray ionization rate 
for H and He (including ionization due to secondary elec- 
trons) in diffuse clouds to be «h£cr ~3x 10~ 17 «h s -1 (see their 
§2.2.2). They also give an approximate expression for the H 
and He X-ray ionization rate £ X r (for the X-ray radiation field 
in the local neighborhood) in terms of x e and the column density 
N w of warm gas surrounding a cold diffuse cloud (see their Ap- 
pendix A). Following Wolfire et al., we adopt Af w = 10 19 cm" 2 . 

For simplicity, we assume that all of the ionization result- 
ing from cosmic rays and X-rays comes from H; 11 in this case, 
x e = x e .o +x H , where x H = «(H + ) /« H - Ionization equilibrium for 
H is expressed by the following equation: 

(&R + 6cR)nH(l-*H) = (a r (B)x c + a g )x H «H • ( 23 ) 
Here we use the "case B" recombination coefficient; i.e., only 
recombinations to states with principal quantum number n > 2 
are included, since recombinations to n = 1 result in an ionizing 
photon which is immediately absorbed by a nearby H atom. We 
take a r (B) = 3.5 x 10- 12 (r/300K)-°- 75 cm 3 s" 1 (Liszt 2001). 



In Figure 8, we display the resulting electron fraction as a 
function of «h for a range of temperatures characteristic of cold 
diffuse clouds (solid curves). We also show the electron frac- 
tion that would result if grain-assisted ion removal were not 
included in the ionization balance (dashed curves). As «h in- 
creases, grain-assisted ion removal becomes more important 
and the results for x e diverge. When n H = 30 cm -3 and T = 
100K, grain-assisted ion removal reduces x e from w 8.0 x 10" 4 
to w 4.5 x 10" 4 , a 45% reduction. 

8. CONCLUSIONS 

The principal results of this paper are as follows: 

1 . We define a rate coefficient a g for the removal of ions 
from the gas phase due to charge exchange with grains, 
including PAHs (eq. 5); "removed" ions either remain in 
the gas phase, in a lower ionization stage, or stick to the 
grain. We compute a g , for several astrophysically im- 
portant ions, as a function of the gas temperature T and 
grain-charging parameter ip and provide a convenient fit- 
ting formula (eq. 8 and Table 2). 

2. We evaluate a critical electron fraction x cr i t ; if x = 
n e /riH < x CI it, then grain-assisted ion removal occurs 
more rapidly than radiative recombination. We find that 
the rates for grain-assisted ion removal and radiative re- 
combintion are comparable in the CNM, but that radia- 
tive recombination dominates in the WNM and WIM 
(Figs. 3 and 4). 

3. We provide formulae for deriving the electron density 
n e from observed neutral/ionized column density ratios 
when grain-assisted ion removal is important, for the 
cases that only two ionization stages are significantly 
populated (eq. 15) and that three stages are populated 
(eqs. 16-21). 

4. We have investigated whether or not grain-assisted ion 
removal can reconcile various estimates obtained by 
Spitzer & Fitzpatrick (1997) for n e in the cold gas clouds 
along the line of sight to HD 215733 (§5.1). Without 
grain-assisted ion removal, the C and Mg determinations 
are consistent, while the Ca determinations are system- 
atically lower. When grain-assisted ion removal is in- 
cluded, the Mg determinations are intermediate between 
the C and Ca determinations, and consistent with both. 
However, the C and Ca determinations remain inconsis- 
tent. 

5. We have investigated whether or not grain-assisted ion 
removal can reconcile the various estimates obtained 
by W99 for n e in the cold gas along the line of sight 
to 23 Ori (§5.2) W99 estimated n e by observing neu- 
tral/ionized column density ratios for several elements 
and assuming that only ionization by the interstellar radi- 
ation field and radiative recombination are important to 
the ionization equilibrium. If grain-assisted ion recombi- 
nation is the only additional process relevant to the ion- 
ization equilibrium, then the various n e determinations 
are not reconciled. However, dissociative recombination 
of CH + may be an important recombination mechanism 

11 The rates for radiative recombination and grain recombination for He + — >He are not very different from those for H + — >H, so the derived electron density is not 
very sensitive to the fraction of the ionizations that ionize H vs. He. A treatment including He would likely result in < 10% increases in n e . 
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for C, and the Na and K neutral/ionized ratios inferred 
by W99 may be incorrect, if their assumed depletions are 
incorrect. Making use of these degrees of freedom, we 
find that most of the observed abundance ratios are con- 
sistent with n e rts 0.03 cm" 3 . However, Mg, Si, and Ca re- 
main discrepant, with inexplicably high neutral/ionized 
ratios. We have considered the possibility that the ob- 
served neutral/ion ratios for Mg, Si, and Ca might be due 
to an additional warm gas component with T w 10 4 K 
and n e > 3 cm -3 . However, such a scenario does not ap- 
pear to be consistent with the observed C + and Si + fine 
structure excitation. Alternatively, the enhanced abun- 
dances of Mg°, Si , and Ca° could perhaps be explained 
if dust is being rapidly destroyed and injecting neutral 
Mg, Si, and Ca into the gas, but it is hard to imagine 
how the dust destruction could occur rapidly enough for 
this to be the case. 



6. We suggest that low temperature radiative recombina- 
tion rates might have been underestimated for Mg, Si, 
and Ca. 

7. We have found that the grain-assisted removal of H + and 
He + can have a significant effect on the electron frac- 
tion in the CNM (§7). For clouds with T w 100K and 
«h ~ 30cm" 3 , x is reduced from 8.0 x 10" 4 to 4.5 x 10" 4 . 
We have assumed that the ionization of C contributes an 
electron fraction of 2 x 10 -4 . 
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APPENDIX 
GRAIN-ION INTERACTION POTENTIAL 

We assume that the electrostatic potential energy of an ion of charge ze at the surface of a grain with charge Ze can be approximated 
by a point charge ze a distance ro from the surface of a conducting sphere with radius a and charge Ze (see eq. [2.9] of Jackson 1962): 

U(Z,z)=^f + 7 ( f e2 ;' - (Al) 
a + r 2ro(a + r y(2a + r ) 

Then, the change in electrostatic interaction energy due to the electron transfer from grain to ion is just AU(Z,z) = U(Z+l,z— 1)- 
f/(Z,z)(seeeq. [7]). 
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Table 1 
Atomic Data 



Species 


IP a 




rv 


r 2 d 


a r+d (T = 100 K) e 


a r+d (T = 6000 K) e 




eV 




s" 1 


s" 1 


3 -1 

cm s 


3 -1 

cm s 


H° 


13.60 


0.37 






8.61E-12 


6.01E-13 


He 


24.60 


0.50 






8.42E-12 


6.71E-13 


C° 


11.26 


0.77 


2.0E-10 


3.3E-10 


8.63E-12 


7.44E-13 


Na° 


5.14 


1.86 


1.3E-11 


1.5E-11 


6.77E-12 


3.00E-13 


Mg° 


7.65 


1.60 


7.9E-11 


7.1E-11 


7.18E-12 


4.66E-13 


Si 


8.15 


1.18 


3.0E-9 


3.1E-9 


9.39E-12 


1.56E-12 


S° 


10.36 


1.03 


7.8E-10 


1.2E-9 


7.46E-12 


5.66E-13 


K° 


4.34 


2.27 


5.6E-11 


6.3E-11 


1.11E-11 


4.16E-13 


Ca° 


6.11 


1.97 


3.7E-10 


4.3E-10 


7.07E-12 


2.40E-12 


Mn° 


7.43 


1.37 


1.4E-10 


1.5E-10 


8.33E-12 


2.16E-13 


Fe° 


7.90 


1.24 


2.0E-10 


2.0E-10 


8.52E-12 


2.48E-13 


Ca + 


11.87 


1.97 


1.4E-12 


2.5E-12 


2.70E-11 


1.02E-12 



a from Lide 2000 
b from Porile 1987 

Ionization rate for the WJ1 radiation field, from Pequignot & Aldrovandi 1986 
ionization rate for the Draine radiation field, from Pequignot & Aldrovandi 1986 

e Rate coefficient for radiative plus dielectronic recombination. a r from rrfit.f (see foot- 
note 6); ay from Aldrovandi & Pequignot 1974 (Na°), Nussbaumer & Storey 1986 (Mg°, 
Si ), Shull & van Steenberg 1982 (Ca°, Ca + ), Arnaud & Raymond 1992 (Fe°). 
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Table 2 



Fitting Parameters for Grain-assisted Ion Removal Rate Coefficient' 



Ion 






Ci 






<-2 




<-3 




(-4 




<-5 


/ ■ 

<-6 




H + 


3.436 


1 


909E- 


6 


1 


.496 


1 


.790E3 


9 


.749E- 


4 


0.5612 


7.245E- 


■5 


He + 


2.175 


2 


.031E- 


■7 


2 


.114 


1 


.134E4 


2 


.197E- 


6 


1.0731 


9.119E- 


9 


C + 


45.58 


6 


.089E- 


■3 


1 


.128 


4 


.331E2 


4 


.845E- 


■2 


0.8120 


1.333E- 


4 


Na + 


2.178 


1 


.732E- 


■7 


2 


.133 


1 


.029E4 


1 


.859E- 


6 


1.0341 


3.223E 


■5 


Mg + 


2.510 


8 


.116E- 


■8 


1 


.864 


6 


.170E4 


2 


.169E- 


6 


0.9605 


7.232E- 


■5 


Si + 


2.166 


5 


.678E- 


8 


1 


.874 


4 


.375E4 


1 


.635E- 


6 


0.8964 


7.538E 


■5 


s + 


3.064 


7 


.769E- 


■5 


1 


.319 


1 


.087E2 


3 


.475E 


1 


0.4790 


4.689E- 


■2 


K + 


1.596 


1 


.907E- 


■7 


2 


.123 


8 


.138E3 


1 


.530E- 


■5 


1.0380 


4.550E- 


■5 


Ca + 


1.636 


8 


.208E- 


9 


2 


.289 


1 


.254E5 


1 


.349E- 


9 


1.1506 


7.204E 


■4 


Mn + 


2.029 


1 


.433E- 


6 


1 


.673 


1 


.403E4 


1 


.865E- 


6 


0.9358 


4.339E 


9 


Fe + 


1.701 


9 


.554E- 


•8 


1 


.851 


5 


.763E4 


4 


.116E- 


■8 


0.9456 


2.198E- 


■5 


Ca ++ 


8.270 


2 


.05 1E- 


4 


1 


.252 


1 


.590E2 


6 


.072E- 


■2 


0.5980 


4.497E- 


■7 



a See eq. (8). 



Table 3 

Inferred Temperatures and Hydrogen Number Densities for Cold Cloud Components Towards HD 215733" 



Component 


T 


»H 




K 


cm" 3 


7 


240 


22 


11 


100 


13 


12 


240 


5 


13 


110 


23 


15 


130 


21 


16 


100 


10 


17 


100 


18 


18 


50 


24 


19 


100 


10 



a From Fitzpatrick & Spitzer 
1997. 



